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The Crab Nebula is the brightest source in the very-high-energy (VHE) gamma-ray sky and one 
of the best studied non-thermal objects. The dominant VHE emission mechanism is believed to 
be inverse Compton scattering of low energy photons on relativistic electrons. While it is unclear 
how the electrons in the nebula are accelerated to energies of 10^® eV, it is general consensus 
that the ultimate source of energy is the Crab pulsar at the center of the nebula. Studying VHE 
gamma-ray emission provides valuable insight into the emission mechanisms and ultimately helps 
to understand the remaining mysteries of the Crab, for example, how the Poynting dominated 
energy flow is converted into a particle dominated flow of energy. We report on the results of 
six years of Crab observations with VERITAS comprising 115 hours of data taken between 2007 
and 2013. VERITAS is an array of four 12-meter imaging air Cherenkov telescopes located 
in southern Arizona. We report on the energy spectrum, light curve, and a study of the VHE 
extension of the Crab Nebula. 
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1. Introduction 

The Crab Nebula and pulsar are two of the best studied astrophysical sources of non-thermal 
radiation in the sky. Their exceptional brightness across the electromagnetic spectrum and relative 
proximity to us make both objects excellent targets to study a cosmic particle-accelerator in great 
detail. 

The central engine of the Crab is a pulsar that provides the energy and particles to power one 
of the brightest pulsar wind nebulae we know. The radio emission like any other pulsed emission 
and particle outflow from a pulsar originates from the vicinity or from within the compact mag¬ 
netosphere that co-rotates with the neutron star and, therefore, appears pulsed. At a distance of 
~ 0.1 pc from the neutron star the pressure of the outflowing cold pulsar wind equals the pressure 
from the nebula and a standing reverse shock is formed|[T], It is in the shock and downstream of 
it where electrons cool via synchrotron radiation that is responsible for the non-thermal emission 
of the nebula from radio up to hundreds of MeV||3|]. 

In this paper we present VERITAS observations of the Crab Nebula over 7 years covering 
observing periods between Fall 2007 and Fall 2014. We present spectral analysis, a search for the 
extension of the nebula, and a study on the variability. 

2. VERITAS 

VERITAS is an array of four imaging atmospheric Cherenkov telescopes designed to detect 
very-high-energy gamma rays between SOGeV and several ten TeV[Q, ^]. The array is located at 
the Fred Fawrence Whipple Observatory in southern Arizona at an elevation of 1280 m above sea 
level. 

The array became fully operational with all four telescopes in Fall 2007. In 2009 one of the 
telescopes had been relocated to its present location giving the array its present diamond shaped 
layout with approximately 100 m separation between the telescopes, yielding an increase in sen¬ 
sitivity of about 30% [^. In Summer 2012 the telescope cameras had been upgraded with high 
efficiency photomultipliers resulting in a sub 100 GeV energy threshold and higher sensitivity. 

3. Observations 

Observations of the Crab are carried out with VERITAS on a regular basis. Besides the primary 
motivation of studying the physics of the nebula and pulsar the observations also provide an easy 
way to monitor the stability of the instrument under the assumption that the Crab is a non-varying 
source. The data are also used to determine the sensitivity of the instrument, verify the Monte Carlo 
simulations. 

Crab observations are carried out with all observing modes of VERITAS, i.e. all zenith angles, 
dark time, moon light observations etc. However, the majority of the data, and this data is used in 
this analysis, are recorded in the Zenith angle range between 8 degrees, where the Crab culminates 
at the VERITAS site, and 35 degrees. 

Data analyzed for this work has been taken between Fall 2007, when all four telescopes be¬ 
came operational, until Summer 2012. In the analysis the data is divided into three epochs, each 


2 



Six years of VERITAS observations of the Crab Nebula 


Kevin Meagher 


representing a significant different telescope configuration, requiring a new Monte Carlo simula¬ 
tion. Epoch one goes from Fall 2007 to Summer 2009, when one of the telescopes was relocated. 
Epoch two spans from Summer 2009 to Summer 2012 when all VERITAS cameras were upgraded 
with high-efficiency PMTs. Epoch fhree sfarfed in Summer 2012 and is sfill ongoing 

The Crab is visible fo VERITAS af Zenifh angles less fhen 35 degrees from Sepfember fo 
March, which are fhe monfhs during which fhe analyzed data had been recorded. All observations 
presented in this paper were taken with a standard wobble offset of O'iS from the source position, 
with the direction of the offset alternating between each of the four cardinal directions. Wobble 
observations allow background to be estimated from off source data in the field of view, canceling 
any sysfemafic effecfs[^, 

For fhe evenf selecfion standard cuts are used fhaf are optimal for such a sfrong source as fhe 
Crab. All resulfs presenfed here had been produced wifh fwo differenf analysis packages in fwo 
independenf analyses. 

4. Spectrum 

The very-high-energy gamma-ray specfrum of fhe Crab Nebula has been measured exfensively 
and measuremenfs are mosfly limifed by sysfemafic uncerfainfies buf af fhe highesf energies above 
lOTeV. The biggest systematic effect in the spectral reconstruction are uncertainties in the energy 
reconstruction that lead to shifts in the reconstructed flux. Af fhe lowesf energies sysfemafic uncer¬ 
fainfies are dominafed by effecfs of fhe frigger elecfronics. 

The reconsfrucfed specfrum is shown in Figure ||. The specfrum fifs well fo a log-parabola 

/ \ —a^p\og(E/E q) 

function =/o ( ^ j - The fit gives a normalization of/o = (3.75zb0.03) x 10“^^ 

TeV^^cm^^s^^ a photon index a = —2.467 ± 0.006, and a curvafure paramefer jS = —0.16± 
0.01. The fit is performed in the energy range of 115 GeV to 42 TeV, the quality of the fit is 
;^^/ndof = 12.9/13. Only statistical errors are considered here, detailed systematic studies are 
ongoing. Extensive studies have been perform on subsets of the data. Only small systematic shifts 
in the spectrum are observed over different cut levels, years, months of the year, wobble offsets, 
and elevation. 


5. Morphology 


The morphology of the Crab Nebula is not yet resolved in gamma rays. It is expected that 
the size of the inverse-Compton nebula is confined fo fhe sfanding reverse shock and fhaf ifs size 
decreases wifh increasing gamma-ray energy O- Af 100 GeV fhe expecfed size of fhe nebula 
is abouf T.5 O- From observafions fhe extension has been consfrained to be less fhan 2'.4 af 
250 GeV p^ ] and less fhan 2'.0 above 1 TeV 112]. 

Wifh VERITAS, like wifh ofher gamma-ray experimenfs, if is nof possible to resolve arcminufe 
sized sfrucfures. Even! by even! angular resolution of VERITAS is approximately 4' af 100 GeV. A 
source wifh an extension of T, resulfs in a subfle broadening of fhe measured spatial disfribufion of 
excess evenfs. In fhe presenf case where a radial exfenf of fhe Crab Nebula of abouf 1'.5 is expecfed 
fhe spatial disfribufion would be 5% broader fhan for a frue poinf source. 
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Figure 1: Left panel; The differential energy spectrum of the Crab Nebula. Right plot: The differential 
energy spectrum multiplied with Data from HESS0 and MAGIC[|^ are shown for comparison. 

Error bars represent statistical error only. 


To perform this analysis, the uncertainty of the point spread function of the telescope must be 
know to better than the expected broadening of the spatial distribution. Monte Carlo simulation 
of the detector proved to be in adequate as it was dominated by systematic errors. Measuring the 
point spread function by observing sources far in the sky from the Crab, such as Markarian 421 and 
Markarian 501, are also inadequate due to the geomagnetic effect. The blazar VER J0521+211[|T^], 
located only 3 degrees away from the Crab was selected as an template to represent the instruments 
point spread function. 

The search for an extension of the Crab Nebula is performed in the parameter 6^, which mea¬ 
sures the radial distance squared between the reconstructed event origin in the sky and the location 
of the gamma-ray source. For the calculation of 6^ it is assumed that the maximum of the spatial 
distribution of excess events coincides with the position of the Crab pulsar in radio: RA = 83.63308, 
Dec = 22.01450 (J2000)[|T^ in case of the Crab Nebula data and that for VER J0521-I-211 the posi¬ 
tion of the maximum coincides with its optical counterpart RGB J0521.8-1-2112: RA = 80.441538, 
Dec = 21.214273 (J2000)|[T|]. 

The 6^-distributions for each data set are fitted with a hyperbolic-secant function by means of 
an unbanned maximum-likelihood method. The radial component of a two dimensional hyperbolic- 
secant function 

= — • sech(V^/w) (5.1) 

KW 

is a good empirical analytical description of the leptokurtic 0^-distribution after background sub¬ 
traction. The parameters of S are and the width of the distribution w. w denotes the 55.1% 
containment radius. 

Because the camera acceptance is not uniform, the background as a function of is not flat. 
When restricting the fit to within 1° of the fit the background can be characterized by a linear 
function. This function is fit by taking background data from the half of the field of view opposite 
from the source and scrabbling it about the center of the camera and fitting the distribution of to 
a first order polynomial. 
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Figure 2: Left: The results of an unbinned maximum likelihood fit of the extent of the Crab Nebula are 
shown in red. The fit was performed using events with reconstructed energies between 160 and 200 GeV. 
Binned data is shown in black for comparison. The Background only hypothesis is shown in green. Lower 
plot shows the residuals. Right: Same plot is shown for the VER J0521+211 which is used as a template to 
measure the point spread function. 


For the actual fit of the excess distribution the following two-dimensional likelihood function 
is defined: 


f{nsig,w) 


-2llog 

02 


^total V ^total J 


(5.2) 


where the summation is over all events, risig, the number of signal events, is used as a nuisance 
parameter, and ritotai is the total number of reconstructed events after event selection. The likelihood 
function is minimized in the range of 0° < 0 < 1° with the MlNUlT software package. 

Figure ^ shows the 6^-distribution of the gamma-ray signal of the Crab Nebula including the 
best fit result averaged over 20 iterations. Figure shows the results of the energy fit in four 
different energy bins. 

Source extension is not evident in the data, so an upper limit on the size of the nebula is deter¬ 
mined with the unified method by Feldman and Cousins [|I^]. For the purpose of constructing con¬ 
fidence intervals, pseudo-experiments are performed to empirically determine the probability dis¬ 
tribution function p{w,s) where s is the source size. For each source size 10^ pseudo-experiments 
or trials are executed. In each trial the origins of the Nsignal events are drawn from a normal distri¬ 
bution in both X and y with a = s. 

The uncertainty in the PSF is included in the limit by determining the PSF for each trial 
from the normal distribution Each event in a trial is given an offset from 

its source location randomly chosen from the radial sech function with the PSF for the pseudo¬ 
experiment. Background events are added by taking events from the Crab data after adding a 
random rotation with respect to the center of the camera. These distributions for the Crab Nebula 
and VER J0521-1-211 are then fit using the same unbinned maximum likelihood fitting code as 
described above. 

Eor the study of the nebula extension the analysis is restricted to post-camera-upgrade data. 
Including more Crab Nebula data in the analysis does not impact the results presented here due 
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Figure 3: The measured width parameter as a function of energy, for both the Crab Nebula and the blazar 
VER J0521H-211. The lowest energy point is used to set a limit on the extent of the Crab Nebula. 


to the dominance of statistical uncertainties in the VER J0521+211 data. Exact values for upper 
limits are still in progress. 


6. Variability 

To investigate the variability of the Crab Nebula we calculated the average flux above 200 GeV 
in bins of one day for every night which VERITAS observed the Crab Nebula. The light curve is 
shown in figure Observations where the elevation was too low to observe the source at an energy 
of 200 GeV were not used in this analysis. No strong evidence of variability is seen in the data. We 
also investigate variability of the flux by comparing the VHE light curve to the one observed in the 
hard x-ray band by Eenn/-GBM[|I^], and in the high-energy gamma-ray band by Fermi-CKV. The 
light curves for these instruments are shown in Eigure Q The correlation between two light curves 
was calculated using the z-transformed discrete correlation function (ZDCE), which is widely used 
to estimate the cross-correlation function of discrete sparse astronomical time series [|T^]. ZDCE 
will calculate correlation coefficients for all time scales up to the duration of the observation. Eigure 
^ shows the results of this analysis. No evidence for correlated variability is seen between any of 
the energy bands on any time scale. The results represented here are consistent with an analysis 
performed by VERITAS searching for correlation specifically during a flare of the Crab Nebula in 
201311 ^]. 
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Figure 4: Light curves of the Crab Nebula as measured by VERITAS, Termi-LAT, and Fermi-GQM.. All 
light curves are normalized to crab units where 1 crab unit is equal to the average of the light curve. 
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Figure 5: Left: Correlation coefficients calculated as a function of a delay between the VERITAS light curve 
and the Fermi-GSbA light curve for the Crab Nebula. Right: Same figure for a correlation search between 
VERITAS and Fenwi-LAT light curves 
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7. Discussion 

In this proceedings, we have shown the results of six years of observations of the Crab Nebula 
with VERITAS. The spectrum is consistent within quoted systematic uncertainties with the spec¬ 
trum published with other instruments. We investigated the morphology of the source and saw 
no evidence of an extension. The lightcurve analysis shows no correlation with variability of the 
source in other wavebands. Analysis of these data including more detailed treatments of systematic 
errors are still ongoing. 
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